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Abstract 

Exact analytic calculations in spin- 1/2 XY chains, show the presence of long- 
time tails in the asymptotic dynamics of spatially inhomogeneous excitations. 
The decay of inhomogeneities, for t — > oo, is given in the form of a power law 



(t/TQ) 



where the relaxation time tq and the exponent vq depend on the wave vector 
Q, characterizing the spatial modulation of the initial excitation. We consider 
several variants of the XY model (dimerized, with staggered magnetic field, 
with bond alternation, and with isotropic and uniform interactions), that 
are grouped into two families, whether the energy spectrum has a gap or 
not. Once the initial condition is given, the non-equilibrium problem for 
the magnetization is solved in closed form, without any other assumption. 
The long-time behavior for t ^ oo can be obtained systematically in a form 
of an asymptotic series through the stationary phase method. We found 
that gapped models show critical behavior with respect to Q, in the sense 
that there exist critical values Qc, where the relaxation time tq diverges and 
the exponent uq changes discontinuously. At those points, a slowing down 
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of the relaxation process is induced, similarly to phenomena occurring near 
phase transitions. Long-lived excitations are identified as incommensurate 
spin density waves that emerge in systems undergoing the Peierls transition. 
In contrast, gapless models do not present the above anomalies as a function 
of the wave vector Q. 
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I. INTRODUCTION 



The spin-1/2 XY chain and its variants are among the most widely used and quoted 
models in theoretical investigations on spin systems. This popularity among theorists is due 
to two facts. On the one hand, this family of models allows for exact theoretical description 
of many static as well as dynamic properties (see Ref. ^ and on the other hand, in some 
important cases, it provides a good description of real systems&H. As an example, we 
note that the dimerized chain, with alternating antiferromagnetic bonds, has been studied 
in relation with the spin- Peierb tra„=itio,fl. and it is tl»ugU to represent the spin degrees 
of freedom of organic compounds that undergo the so called Peierls distortion. The present 
state of the art in fabricating low dimensional systems, with the material science technology 
developed after the synthesis of the superconducting cuprate oxides, may now tailor com- 
pounds that reveal a wealth of new magnetic phenomena, including random spin chains, 
spin ladders, and doped magnetic chains among other systems with exotic behaviors. For 
theorists, these systems are fascinating, with no parallel in classical or three-dimensional 
physics. Being systems of low dimensionality with low values of the spin, they are domi- 
nated by quantum effects^. Compounds such as PrC/3, PrEtSO/^ {Pr Ethyl Sulfate), and 
CS2C0CI4 are among the quasi-one-dimensional systems, whose low temperature properties 
are thought to be described by the XY model. 

A great deal of theoretical information has been gathered concerning equilibrium prop- 
erties, with calculations of quantities such as the specific heat, the magnetic susceptibility, 
and equal-time spin-spin correlation functionsEI'ii'El. 

As for the dynamic properties, the quantities most thoroughly studied are the time- 
dependent spin-spin correlation functions (TDCF) 

Sf{t)St) = (exp (^m) ^;exp (-^m) St) , (1) 

where H is the Hamiltonian of the system. is an index for the spin component (/i = x, y, z), 
and < ... > is the equilibrium averageSBBil. They are important for the description of such 
dynamic phenomena as magnetic resonance, magnetic neutron scattering, spin diffusion, 
and other relaxation propertieflil. It should be noted however, that their application 
is restricted to situations near the equilibrium state, where the linear response theory is 
validiiil. In particular, for the uniform XY model, an old calculation for the time dependent 
autocorrelation < SQ{t)SQ{0) > showed the absence of spin diffusion in the limit t oci. 
This behavior was thought to be accidental and specific of the XY model for spin 1/2 (see, 
for example, Ref. However, recent research has shown that this surprising property 

is shared by a whole family of integrable models, and is attributed to the existence of a 
macroscopic number of conservation lawS. This remarl^ble result has been probed in 
magnetic resonance experiments in some ID systemsH. We will turn to this point later. 

In spite of the wide range of their applicability, we note that TDCF do not give a direct 
description of the time evolution of the non-equilibrium process00. This description is 
achieved through a different method, that we briefly discuss below. 

This second direction in investigating the dynamics relies on calculation of averages of 
the type: 

{A\ = Tr [p{t)A] , (2) 
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where p{t) is the density matrix which satisfies Liouville equation 

th^^p{t) = [H,p{t)] . (3) 

Averages of such a kind give a direct account of the non-equihbrium evolution of the physical 
observable A, independently of how far the initial state is from the equilibrium state or from a 
stationary one. Unfortunately, calculations of these quantities are much more involved than 
the calculation of TDCF (|l]) for linear response, and few exact results are knownlli'i&^l. 

Whether the dynamical process is near or far from equilibrium, most works in the liter- 
ature deal with cases where the initial state is spatially uniform. This premise is assumed 
in explicit or non-explicit way, and the methods developed to solve the problem are heavily 
based on it. In contrast, the study of spatially nonuniform excitations is practically new, 
in spite of its interest, both theoretical and experimental. This problem is important for 
a deeper understanding of dynamical processes in many-particle spin systems with strong 
exchange interactions. Inhomogeneous initial states can be prepared in real systems by ex- 
ternal actions, for instance, strong inhomogeneous magnetic fields or acoustic waves. On 
the other hand, from the theoretical side, exact results on the dynamics (as in the case of 
the XY model), can elucidate details of spin-spin relaxation processes in more complicate 
systems. 

In the present contribution, we will adopt the method based on formulae @ and and 
will analyze in detail the long time evolution of the magnetization in three versions of the XY 
model. The initial excitation is always prepared in the form of a Spatially Inhomogeneous 
Magnetization (SIM), and the calculation is done in exact analytic form. The three variants 
of the XY model that we consider are enumerated below: 

I Isotropic dimerized XY model with Hamiltonian: 

N M 

Hi = —psgh^ Sj — ^ [Jl{S2m-1^2m + Slm-lSlm) + 
j=X m=l 

(4) 

, J ( qx qx j_ qy qy \i. 



II Isotropic XY model in a staggered magnetic field with Hamiltonian: 

M N 
m=l j=l 



III XY model with bond alternation without magnetic field: 

M 

-^3 = ~ X] (-^1*^2™- 1*^2™ + -^2'S'2m-l'S'2m + -^2'S'2m'S'2m+l + -^I'S'lm'S'lm+l ) • (6) 
m=l 
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We will refer to them as models I, II, III, respectively. In expressions (^-®, N = 2M, 
/is is the Bohr magneton, g is the gyromagnetic ratio, h, hi, h2 are magnetic fields, J, Ji, 
J2, /i, I2 are exchange integrals, (Ji > J2 > 0, and J > 0), and cyclic boundary conditions 
are assumed. The above restriction on Ji, J2 and J is not important in our problem and is 
introduced for convenience. 

To investigate SIM dynamics in the above models, we will use a previously developed 
method given in Ref. This method was applied earlier by one of the authorsBB, to the 
same study for the isotropic XY chain in a homogeneous magnetic field, with Hamiltonian: 

N N 

We will refer to this model as model IV subsequently. We remark that highly nontrivial 



results were obtained in Ref. |3^ , |35| . In particular, it was shown that in the t ^ 00 limit, some 
of the spatially inhomogeneous excitations do not disappear and are still time-dependent. 
The time evolution of SIM's (and also their spatial distribution) is probed through the 
computation of the Fourier components (^Sqlt)'^ of the magnetization as a function of the 
Q-wave number. The calculation yields a relaxation process in the form of a power law 

where the exponent u depends on the initial state. In Ref. ^ it was also studied the 
anisotropic XY model in the limit of strong anisotropy. It was found that the exponent 
u in changes discontinuously at some critical values Qc of the wave vector Q. The 
values Qc depend solely on the parameters of the Hamiltonian and are not connected with 
the preparation of the initial state, nor the particular component of the magnetization 
that is relaxing. Moreover, in the limit Q Qc, the relaxation time r of (H) diverges 
to r ^ 00, with the corresponding slowing down of the process. This phenomenon is at 
variance with some conventional views for the time evolution of physical quantities in many 
particle systems, according to which all quantities must be temporal independent in the 
limit t ^ 00 . The concept of spin temperature, widely used in non-equilibrium magnetic 
phenomena@, is based on the assumption that the spin-spin relaxation is much faster than 
the final relaxation to the lattice degrees of freedom. 

Model I was also preliminary investigated by the some of the authors, and similar con- 
clusions were attained for the dimerized XY modelB, in the sense that a critical-like slowing 
down of the relaxation process takes place at special points ±Qc of the Q-space. In other 
words, for | Q |— > Qc, the inverse time scale for some of the oscillating components of 
SIM, goes to zero following the power law ~ || Q | —Qd- Such behavior is not surprising 
for the case Qc = 0, for which the corresponding value of the total magnetization is a 
constant of motion, but is unusual for Qc 7^ 0, where the corresponding Fourier component 
Sq is not conserved. 

We note that Q = is the only critical point for the uniform isotropic XY model (model 
IV and limit Ji = J2 for model I and hi = h2 for model II ). The t —>■ 00 behavior is 
dominated by one oscillating component, and this component has no critical properties. 
The present paper is devoted to elucidate this remarkable difference within a more general 
context, by the extensive study of a whole family of models. We suggest that the dissimilar 
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behaviors are due to the presence of a gap in the spectrum of the energy excitations, the 
uniform isotropic model being gapless. 

Our paper is organized as follows: In Section II, we outline the main steps in the analytic 
calculation. All the models are diagonalized by means of a modified Jordan- Wigner trans- 
formation, which maps the spin model into an equivalent fermion Hamiltonian. Then, the 
average < Sq >t is calculated in exact closed form, and its asymptotic behavior for t ^ oo 
is obtained after taking the thermodynamic limit —>■ oo. In Section III, we give a detailed 
analysis of the results, including all models and the long time behavior of the magnetization. 
Section IV closes the paper with the final discussions. 



II. MAIN STEPS OF SOLUTION 

The quantity of our interest is the time-dependent z-component of SIM 



) = Tr 



Pit)S^Q\ , (9) 
where p{t) is the density matrix of the system, and 

N 

^Q = E S^eM^Qj) , (10) 

is the Fourier transform of the magnetization. In (P), Q is the wave vector characterizing 
the spatial inhomogeneity of the initial state {Q = 27m/N, n = —N/2 + 1, N/2). To 
calculate < Sq >t one can use the identity: 



where p(0) is initial density matrix, and Sqlt) is the spin operator in the Heisenberg repre- 
sentation 

5^(t) = expQm)s^exp ^-^Ht 

To obtain exact results for < Sq >t, for each of the models considered, we follow the steps: 

i) Diagonalization of the corresponding Hamiltonian using the fermion representation of 
the Jordan- Wigner transformation. At the end, we get a free-fermion system; 

ii) Transformation of 5*^ in terms of fermionic operators to get the time dependent ^'^(t); 

iii) Averaging of SQ{t) with the initial density matrix p(0) . 

The last step depends on the form of initial matrix p(0). As it was shown in Ref. |3^, 
exact solutions of this problem may be obtained in the case when p(0) is a functional of only 
one component of the spin operator: 

p(0) = F(^^), (7 = x,y,z). (11) 

Let us note that such initial state can actually be prepared in real systems at low tem- 
perature, with a strong nonhomogeneous magnetic field directed along the coordinate axis 
7- 
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A. Diagonalization 

Methods for diagonalization of Hamiltonians (iD-(|D are well known (see for example Ref. 
P, p!T| , p!3D ■ Minor differences are specific for the version of the model to be solved. In our case 



we have employed the following procedure: 

a) Use of the Jordan- Wigner transformation to change from spin to Fermi operators 

S^ = LM + bj)/2, S] = LM-b,)/2t, S; = b]b,- 1/2, (12) 



where 



= - 1), Li^l, (^2 = 1); 



1=1 



b) Introduction of two types of Fermi operators for even and odd sites: 

C = &L, C = &L-i, m = l,2,...,M, a = ±l, (13) 



and of their Fourier transforms , ci^ : 

M M 



exp{—iakm)c'^, d'^ = —= ^ exp[—iak{m — l/2)](i^, (14) 



k = 27rmi/M, with mi = —M/2 + 1, ...,M/2 , where we have adopted the compact 
notation 



c = c , = , 

d~ = d , d'^ = d^ ; 



(15) 



c) After the above steps, Hamiltonians (^-(§|) are transformed into quadratic forms in 
terms of operators (c^, c?^ ). Diagonalization of these quadratic Hamiltonians is stan- 
dard, but depends on the specific model under consideration. For Hamiltonians (^ 
and (^ it can be done with the help of the Bogolyubov transformation given below: 

ct = uuVt + vWt , dl = v^'^Vk - uuPt , (16) 

where(r7^, /3^) are new Fermi operators, f^" = (f^)*, and Uk is a real function of k . 
We are using the same convention ([T5|) for a = ±. 

For Hamiltonian (^), the functions (ufc, Vk) are given by: 

Mfc = , < = exp(ia^fc) (17) 

where 
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1 

TT6 



tan {k/2) , 6 = J2/J1 < 1 



For Hamiltonian (|^), one gets: 

Uk = [{Pk - cUrf)/2Pfc]i/2 , Vk = [{Pk + uJd)/2PkY/hign{cosk/2) 

where is real (independent of a) and 

uJd = fJ^Bgih - h2)/2h , Pk = {ujj + uj^cos'^ k/2Y^'^ , Ue = J/fl 



d) Following Ref. |Tl], diagonalization of Hamiltonian is achieved through: 



(18) 



(19) 



(20) 



As a final result, all the three Hamiltonians can be represented in the diagonal form 



Hi = fiT,k 
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(z = 1,2,3) 



(21) 



except for constant terms that are not important for the dynamics. The index i = 1,2,3 
refers to models I, II and III, respectively. The dispersion relations of (^) are given by: 



(1) 



-ujo + uJiRk , ojf^ = -uJs + Pk , i^k' = -{h/K) cos {k/2) 
-uo - uJiRk , el') = -Us - Pk , = ih/h) cos ik/2) 



,(3) 



e 



(1) 



where 



Rk = VI + 6^ + 26 cos k , Pk = Jul + ul cos^ {k/2) 



= fj'Bgh/h , Us = fiBg{hi + h2)/2h , ui = Ji/2h . 

Note that here and in the following, we neglect, as usually done^'ii'El, the boundary term 
of the order 1/A^ which appeared in a chain with cyclic boundary conditions. The thermo- 
dynamic limit N ^ 00 will be taken at the end of the calculation. 



B. Calculation of < Sq >t 

Using formulae (ID) and (|r^)-(n), one can easily obtain the following expression 



k 



(22) 



where p = k — Q, q = k + Q, with definitions common for all models. To obtain the time 



dependence of Sq, we first express the relation (|2^ ) in terms of the canonical r/'s and /5's 
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through transformations ([Tq) for models I and II, and (^) for model III. Then, we insert 
the time evolution of such operators 

Vtit) = Vk exp(zaa;f t), P^{t) = exp(zaG«t) , 2 = 1,2 

to get 

S^Qit) = -M5q,o + 

+ Ek {{upUq + v*Vq) exp[i{Qi{k, Q)t]r]^T]q + i^pUq + VpV*) exp[-iQi{k, Q)t](3^(3q+ 

+ {UpV* - UqV*) exp[inf{k, Q)t\r]^ f3q + {UqVp - UpVq) exp[-2l]+(/c, Q)t](3^T]q^ (23) 

for models I and II, where 

n+(fc,Q) = ^«-e», ^{k,Q) = uf -ujf, 2 = 1,2 

and functions ( Up-, ^p) ci-rG defined by formulae (|T7p and (|I^) respectively. 
In the same way, for model III, we obtain: 

Sm = E {(^-X exp[2(a;f ) + Qf)t] + r^_A exp[-2(^f ^ + ))t]} . (24) 

k 

To calculate the time dependent average < Sq >t, we use the identity 

which means that, equivalently, one can calculate the average of operator Sqlt) in the 
Heisenberg picture in relation to the initial density matrix p(0). 

From equations (|23|) and (plD, one sees that this problem reduces to calculation of av- 
erages of the type < rjj^riq >o , < P^Pq >o including all the combinations of canonical 
operators. The calculation is described in detail in Appendix 0. 

Below, we will summarize these results together with the ones obtained in Ref. 35 for 
model IV. We remark that exact solutions can be obtained for special forms of the initial 
density matrix p(0). We assume here that p(0) is a functional of only one component of the 
spin operators, i.e., p(0) = F{S^), (7 = x,y,z). 

1. In the case 'j = x (the case 7 = ?/ is the same by symmetry) we have 

{Sq)^ = (25) 

for models I, III, IV, and 

■ 5^ X = -3c.ee., sin(Q/2)e^«/^ ( >o 

[cos[m^ {k, Q)] - cos[m^ {k, Q)] } dk, (26) 



— I 
TT 

^ sin^(A;/2) 



P P 



for model II, where 



N 

(4)o = E (^J5jVi)^exp(2A:j) (27) 
is the Fourier transform of the short-range order correlation of the x component. 
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2. In the case 7 = we have different expressions for all the cases, 
(a) For model I, we get the formula: 

Sq)^ = ^ f {F+{k,Q)co^[tn+{k,Q)]+F-{k,Q)cos[m^{k,Q)]]dk, (28) 



with 

F^ik.Q) = A^{k,Q) < S5 >„ +B^{k.Q) < >„ , 

(b) For model II, we have the more explicit form: 

{S'Q\ = i^{Sc^)j-^dk 

{-cj2 cos(p/2)cos(g/2) /PpPg cos[tn^(fc, Q)]+ (30) 
+ [1 + {ujj + ul cos(j9/2) cos(g/2)) /PpP,] cos[t^]2 (A:, Q)] } • 

(c) For model III, the compact expression: 

sh), = ^{sh >o £ Q)\dk (31) 



where 



VL{k, Q) = vq cos(fc/2 + aq) ; 

tanag = — ^ tan(Q/2) , (32) 
h — 12 



hvQ = ^Jlf + li - 2/1/2 cos Q . 

Let us note that at the specific value Q = n, the time dependence of (•S'g^^ is 
given by the simple formula: 



( SJ ), = ( S; )„ J„ ( 



h 

where Joix) is the Bessel function of first kind. 



(d) For model IV, it was found in Ref. |35 



Sh)=(Sh)^M^Qt), (33) 



with rig = 2{J/h) sinQ/2 . 
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III. ANALYSIS OF RESULTS 



In this section, we firstly discuss general properties of the time evolution of the magne- 
tization, which can be inferred directly from equations (|A5|) and (p5D-(^). In addition, we 
will study in detail the asymptotic behavior at long times, in the limit t ^ oo. We will show 
that the evolution of SIM's displays interesting unusual features, depending on the wave 
vector Q characterizing the inhomogeneous initial state. 



A. General properties of < Sq >t 

Let us first consider model I in the case 7 = 2; for the initial density matrix, since it 
presents a specific feature that is absent in the other models studied here. From equations 
(^) and (p^)-p9|), one sees that the value of the magnetization < Sq >t at t > 0, depends 



on the initial values of the two Fourier components < Sq >o and < Sq^^ >q. This means 
that, if only one component of Sq exists at t = 0, for example for wave vector Qq, i.e. 

< Sq >0 = < Sq^ >o 5q,Qo ) 

then a twofold response develops for t > 0, with components < Sq^ >t and < S^_^_q^ >t, 
whose time evolution is described by equation (|A5|) , with F^{k,Q) = A'^{k,Qo) < Sq^ >o 
5q,Qo and F^{k, Q) = B^{k, Qq + ti) < Sq^ >o ^q^Qq+tt , respectively. In contrast, for models 
II-IV, if there is only one initial component < Sq^ >o, only one component < Sq >t of the 
same wave vector Q = Qo exists at subsequent times, t > . 

When the initial condition is prepared with 7 = 0;, the average < 5*^ >t vanishes 
identically for models I, III and IV, for any t > 0. This result may be understood intuitively, 
without any calculation, if one notes that the part of the Hamiltonian that describes the 
coupling of spins with the magnetic field, commutes with the other one describing the 
exchange interactions (in the case of model III, this statement is trivial for we only have the 
exchange part). Denoting the latter as H^x, we get the time evolution 



S'n) =Tr 



p(0) exp QiJex. S'q exp (-^^e. , (34) 



which is governed by the exchange only. Inside the trace in (3|), we now perform a unitary 
transformation consisting in a rotation by vr around the x axis. Under this operation, the 
spin operators transform as follow: 

QX QX oy cy QZ QZ 

, Dj- Dj- , . 

The exchange Hamiltonian H^x is invariant under this transformation because it is a sum of 
products of the form SjSj_^-^ and SjSj.^^. The initial density matrix p(0) is also invariant 
since it only depends on the Sx component of the spins (7 = x). So, the right hand side of 
equation (|3^) will change its sign, yielding < Sq >t= — < Sq >t i.e. < Sq >t= . 

In the case of model II, the exchange part does not commute with the total Hamiltonian 
This leads, in general, to a nonvanishing < 5*^ >t, for t > 0, which results proportional 
to < £q^t, >o 5 where < el >o is the Fourier transform of the spin-spin correlation function 
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< SjSj_^^i >o evaluated at the initial condition (see relation (]27|)). If the initial state is 
prepared with a single Q-vector in the form < ef. >o=< >o ^k,Qo , then < •S'^+Qq >t is 
the only component that exists for t > . 

For Q = Qo = (homogeneous initial state), and for any 7, we get the (evident) result: 

< Sq^q >t = < Sq^q >o 

for models I, II, IV. This is a direct consequence of the conservation of the total 2;-component 
of the magnetization. In the case of model III, when 7 = 2;, we get the closed result: 

< Sq^q >t=< Sq^q >o Joi^ot) 

where h = \Ii — ■ 

Another important property valid for models I, II and IV, and for any 7, is that the 
evolution of SIM's is independent of the value of external field h . This fact is not surprising 
because the Zeeman term of the respective Hamiltonians commutes with both, the remaining 
part describing the exchange interaction and the operator 5*^ . This way, the time evolution 
of SIM's is determined by the exchange term only. 

Special important limits are 5 = 1 (Ji = J2) in model I and hi = h2 in model II, that 
map onto model IV. In both of these cases, the corresponding formulae for < 5*^ >t reduce 



to (BSD which was obtained in Ref. |35. 



Unfortunately, it is difficult to get more information from formulae ( [A 51) , (|26|), (p9[) , and 



(PPP, because they are rather involved. So, in the following, we will consider the limiting 
case t ^ 00, that will be studied using the stationary phase method. 



B. Long time behavior of < Sq >t for the initial condition 7 = 2 

1. Model I 

According to the stationary phase method^, the long-time evolution of < Sq >t is dom- 
inated by the contribution of the stationary points of the functions ^lf{k, Q) . The number 
of these points depends on the value of the wave vector Q. For convenience of further dis- 
cussion, let us introduce here the so called critical values of Q : 

Qd = arccos5, Qc^ = vr — arccos5, = 0. (35) 

The above values are determined in a standard manner through the stationary phase method, 
being the locus where stationary points become degeneratec3. Let us first examine the role of 
Qf{k, Q). It is easy to show that it has five nondegenerate stationary points in the interval 

Qci "^l Q \^ Qc2 • 

h = 0, k2 = -fcs = TT, 

k4 = -fcs = (fiQ, (fiQ = arccos(-5"^ cosQ). 
In turn, for Qc^ <\ Q |< vr and Qc-^ <| Q |< Qc^, one gets three nondegenerate points at: 

ki = 0, /C2 = -k'i = IT. 
Exactly at the critical values, we obtain: 
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i) For Q = ±Qci, one nondegenerate point ki = and two degenerate ones ^2 = —ks = n 
{pt{±TT,Q)]" = [nti±7i,Q)]"' = 0, [fi|(±7r,g)](^^) ^ , where derivations are taken 
with respect to k); 

ii) For Q = ±Qc2 , two nondegenerate points ^2 = — fcs = vr, and one degenerate at 

k = k^ = o i[ntio,Q)]" = [ntio,Q)r = 0, [fi|(0,Q)]^^ ^ O); 

iii) For Q = Qc.^ = , there is no time evolution, since the magnetization < Sq^q >t is a 
constant of the motion, 



The function {k, Q) vanishes at Q = Qcs = and Q = ivr. For any other wave 
number Q, it has two nondegenerate stationary points: 

k^ = —k-j = ipQ, 'ipQ = arccos(— (5cosQ). 

This finishes the analysis of the stationary points that dominate the long-time behavior of 
(p8|). The corresponding asymptotic development of < Sq >t for t — > oo, can be represented 
as a sum of several oscillating components with Q-dependent frequencies and amplitudes: 



S'Q)^r^j:S.iQ,t), (36) 



S,iQ,t) = J2a^,iiQ) - cos[tn,{Q) + a,,i] , (37) 



=0 



where Tj = Ti{Q) are functions of Q and the exponent z/j assumes the values 1/2 or 1/4 (the 
latter value will be discussed in detail below). The number of components Si{Q,t) depends 
on Q. For Qci <| Q \< Qc2 (excluding Q = ±7r/2, which is a degenerate point where Qi and 
^^2 coincide), there are four components. For <| Q |< Qci and Qc.-, <\ Q \< tc (excluding 
accidental degeneracies), there are only three terms in (^). The frequencies ^li{Q) and the 
inverse of the characteristic times r^"^ of (|37D are given below: 

ni{Q) =tui {1 + 5^ + 25cosQ 



1/2 



T^=uJid{l + dcosQ) -3 , Q ^ ±Qc 

n^iQ) = cji (1 + (5^ - 25 cos g)^^^ , 



r2 ^ = UJ16 {1 - 6 COS Q) I ^ ^'^^ ^ I ^ Q ^ j^Q 



(38) 
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^3 ^ = ^loTT^TaT (^'^Q-Q + R^q+q) I sin Q I sin^'g, Q 7^ Qcg = 



^iiQ) = uji 



2(1-6 , 



2 ; ' 



C2 



We display typical Q-dependences of fli{Q) and r^'^in Figs.l and 2, respectively, for partic- 
ular values of 5. Slowing down of the relaxation process occurs at points where the ^""^'s 
vanish. 

At a generic Q-point, the first nonvanishing terms in (|37|) have the form: 



t \ -1/2 

SiiQ, t) ^ ai,o(g) ( - cos [tn,{Q) + a,,o] , = 1, 4), (39) 



with the amplitudes 

a,o(Q) = al,{Q) {S^q)^ + zal,{Q) { 5^+^ )^ , 

showing that the value at Q depends on the two initial components (^S'g^^and (^SQ_^_Jj^. 
Formulae for the real and imaginary parts a^Q^Q) and a^Q^Q), as well as the phases ai^, are 
given in Appendix 0. 

The dynamical process described by (|36D and (^), with the explicit formulae given in 
(PSP and Appendix 0, is remarkable, since it exhibits long time tails in the relaxation of 
SIM's to the spatially homogeneous state. In contrast to the exponential relaxation, which 
is characterized by a single parameter that yields the time scale or the relaxation rate, the 
power law relaxation given by (|37D and (pOj) is characterized by two parameters: rf^, which 
determines the inverse time scale of the process, and the exponent z/j, which determines the 
relaxation rate. In general, no conservation laws or long-lived hydrodynamic modes seem to 
be associated with the above long time tails. 

In the neighborhood of critical points, the relaxation of some of the components of 
begins to stop, with the corresponding relaxation time diverging, as shown in Fig. 2. 
Exactly in the limit, the corresponding exponent Ui jumps discontinuously from 1/2 to 1/4. 
This slowing down of the relaxation process is very similar to the critical slowing down found 
in phase transition phenomena.0 

Let us summarize below the singular behavior of the relaxation time at special points 
(they are displayed in Fig. 2): 

a) rf ""^ — i> when Q ±<5c2; for any S ; 

b) when Q iQci, for any 5; 

c) T^^ when Q ±<5ci, ±Qc2) foi' any 5 ; 

d) T^^ when Q ^ 0, ±7r, for any S. 
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All the Tj ^ vanish according with the law: 

~ 1(1 g I -Qc)\ for \Q\^Qc. 

We now give the behavior of components of (|36D at the critical points: 

I In the limit Q ±<5ci, the components S^iQ^t) and S^iQ^t) merge into a single 
component with v = 1/4: 

52(±Qci,t) = a2,o(±QcJ (:^) '^'cos[tn2(gcJ + vr/8], (40) 

with 

a2,o(±gcJ = |?r(l/4) [(1 + 6) {SIq^X + ^Vl^'{Sl^Q^X 
r-i = 3a;i5V(l-'5'), 
where r(x) is the Gamma function. The other components follows the law (pQ]) with 

II Analogously, for Q ±<5c2) the components Si{Q,t) and S'4(Q,t) join together with 
= 1/4, and the time dependence: 

/ t \ -^1'^ 

^i(±g,,,t) = ai,o(±gcJ - cos[^^]l(g,J + 7r/8] , 



where 

«i,o(±gcj = |pr(i/4) [(1-5) (^iQ^^)^+^v/r^(^^iQ^j^; 

The other components follows the law (^) with Q = ±Qc2] 

III The component 5*3 (Q, t) is critical at Q = Qc-j = 0, with the relaxation time ~ 1/ |<5| 
diverging in the vecinity of Q = 0. Note that the frequency vanishes in the limit 
Q —>■ 0. The criticality of Qc-j has a different connotation, since it is related to the 
conservation of the total 2;— component of the spin 

N 

QS \ ^ QZ 

'-'q=o — 2-^ ' 

which is valid for all the model with axial symmetry along the axis. The divergence 
of r3 near Q ~ is connected to the stability of hydrodynamic excitations of long 
wavelength; 

IV Note the exception of Q ^ ivr for the component 6*3 (Q,^), because the amplitude 
'23,o(Q) goes to zero in this limit. 
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2. Model U 



When applying the same method to get the asymptotic behavior of ( |30D for long times, 
one sees that the work is facilitated by the development of the previous section. In fact, the 
frequencies f2f (fc, Q) and the functions that appeared in ( pOD of model II, can be mapped 
onto the corresponding fif and uJ\Rk of model I, with the substitution: 



«2,o(g) = ( SI 
«3,o(g) = ( SI 

«4,0(Q) = ( SI 



" ■ / o , , , ' ^1 ^V^d + ^e + ^- (41) 

\/^d + ^e + \^d\ ^ ^ 

To avoid numerous definitions, we use the previous symbols b and io\ throughout this section, 
with the meaning given above in (0). With these variables, the stationary and critical 
points, the inverse characteristic times and the frequencies coincide with those of 

model I. 

The expression for {Sq^^ has the form (|37|) with the following values of the first nonva- 
nishing amplitudes ai k{Q) (at a generic Q-point): 

5(1 -5)2(1 -cosQ) 

^ 5(1 - cos Q) 

- Rl^. ' 
2 1 + 5 cos Q 
n T+~6 ' 
2 5 + cos Q 
T+5~ ' 

The corresponding phases ai^i are given by formulae (^^ ). Note that the component with 
frequency ^i{Q) decays to zero as t"^/^ (exponent u = 1/2 and / = 1), i.e. more rapidly 
if compared with the similar of model I . We now summarize the behavior at the critical 
points: 

i) At the critical point Q = ±Qci, the components S2{Q,t) and S4^{Q,t) are degenerate 
and decay to zero with exponent u = 1/4, following the law (^OD with 

\/24 5 
a2,o(±gci) = Vr(l/4)- 



27r ' ' ' 1 + 5 ' 

ii) At Q = ±Qc2, the components Si{Q,t) and Si{Q,t) merge into a single component 
with the evolution (z/ = 1/4 and / = 1 in (HO)) 



S^{±Q,„t) = a^,^i±Q,,) {-) cos[tn^{Q,,) + n/8] (42) 



with 



ai,i(±Q,J = ^{3/2f/' r(3/4) 



27r' ' ' ' ' '1 + 5 
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3. Model in 



The long time behavior of ( ^I]) is much simpler than those displayed by models I and II, 
for the phase function VL{k, Q) in (^) has only one stationary point at /c = /ci = — 2q;q . In 
the asymptotic limit t — > oo , it may be represented in the form 

where Vq is given by formula 



32) and r, 



Q 



C. Long time behavior of ( Sq ) for the initial condition 7 



As it was mentioned in Sec. [II B| , \Sq^^ , for the initial condition 7 = x, vanishes 
identically for models I, III, IV, and only has a nonzero value for model II. For the latter, 
the asymptotic expression of < 5*^ >t for t —>■ 00, has the form of (|37|) with the quantities 



fij and Ti defined through formulae ( p8| ) and (|iT|). The coefficients ai^i{Q) for the first 



nonvanishing amplitudes in (|37|) , at generic Q-points, are: 



V27r 

a3,o(g) = - ( )^ ^^^^^^ (1 - e^^) (1 + 5 cos g) , 

«4,o(g) = ( 71=^7^^^ - "^''^^^ + ' 

where ( )q is the correlation function defined in (|27|) . We remark that the first nonvan- 
ishing term for the component with frequency ili has u = 1/2 and I = 1, yielding a t"^/^ 
power law for the time decay. The corresponding phases are given in ( P3|) in Appendix p. 
We make a summary of the critical behavior at singular points: 

i) At the critical point Q = ±Qci, the components S2{Q,t) and S4{Q,t) merge into 
one component with time evolution given by (^OD, with u = 1/A and / = 0, with the 
amplitude 

«2.o(±gcj = |pr(i/4)^(i-5T^v/r^) ; 

ii) In analogous form, at Q = ±Qc2 the components Si{Q,t) and S4,{Q,t) become degen- 
erate with the time evolution of (0), with u = 1/4 and I = 1, and the amplitude 

aiA±Qc.) = ^(3/2)'/' r(3/4)v^ ^1 t ^/f^j • 
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IV. FINAL DISCUSSION 



We have studied the relaxation to the homogeneous state, of an initial excitation which 
has been prepared with a nonhomogeneous magnetization profile along the magnetic chain 
(SIM). The time evolution of this excitation is probed through the calculation of the Fourier 
component (^SQ{t)'^ of the magnetization, which is analyzed as a function of the wave vector 
Q. We use periodic boundary conditions for the spin Hamiltonian and take the thermody- 
namic limit (A^ oo) before studing asymptotic long times. As remarked in Ref. ^, the 
order of the limits is very important. All the models treated here, show long time tails in the 
relaxation of (^^^(t)^, which is apparent from the asymptotic study for t ^ oo, developed 
in the previous sections. This behavior manifests itself in the form of a power law decay 
in the long time evolution of SIM's. This slow relaxation is a remarkable result for its own 
sake and is probably due to the absence of dissipation in the models. Since the systems are 
isolated at zero temperature, the dynamics is exclusively driven by quantum fluctuations. 

Our calculation also shows striking differences, when one compares the behaviors of 
models I and II on the one hand, with models III and IV, on the other. The dissimilarity 
lies in the presence of critical values for the wave vector Q , where we get a slowing 
down of the relaxation process. Near critical points, the time scale Tj , giving the relaxation 
of a specific oscillating component of (^^^(t)^, diverges with the law ~ || Q | —Qal-, for 
I <5 Qc,, reducing the damping of these excitations. Such behavior is not surprising for 
the case Qcs = 0, for which the corresponding value Sq^q is an integral of the motion, but is 
unusual for the case Qc^ ^ 0, when Sq is not conserved. In this sense, the criticality at iQci, 
±Qc2 has different implications in the theory. The position of the critical values iQci and 
±Qc2 depends only on parameters of the Hamiltonian, and in the limit Q —>■ iQd ? =tQc2 ? the 
exponent u jumps discontinuously, the relaxation rate becoming slower than for excitations 
with \Q\ 7^ Qci2- This means that the critical components will be the only surviving ones 
at long enough times. These long-lived excitations may be pictured as spin density waves 
(charge density waves for the fermion model) with incommensurate wave vector Qci 2 ) which 
are developed in our gapped models as a consequence of the Peierls instability^. In contrast, 
the dynamics of models III and IV does not have such anomalies. There is no gap in the 
spectrum, and in the language of particles, the systems are always metallic and display 
no critical points (except the point Q = 0). In the asymptotic regime t —>■ (yo , only one 
oscillating component exists which displays no critical behavior. 

Peculiar properties of the energy spectrum seem to determine the critical relaxation 
phenomena. In fact, models I and II display an energy gap between the ground and excited 
states, whose size is monitored by the parameter 6 (for model II through the transformation 
(PI)). The existence of critical points is directly related to the presence of the gap, since the 
latter changes the curvature of the dispersion relation for the energy. If one looks at formula 
(P3|) for the Fourier component Sglt) of the magnetization, one realizes that processes which 
contribute with the ifij^ frequency (the one that is critical) come from transitions between 
both branches of the spectrum (interband transitions), where a particle with momentum 
q = k + Q is destroyed in a given band, and a particle with momentum p = k — Q is 
created in the other one, with a net momentum transfer Ak = ±2 \Q\. So the criticality 
at Qci,2 is not directly related to the one particle spectrum, but to a big density of states 
for those interband transitions that exchange the same momentum Ak = ±2 \Q\ coherently 
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(constructive interference). This nesting condition is attained at Q = Qci2) with the onset 
of a spin density wave for the magnetization. 

However, we have found that the above anomalous behavior in the relaxation is not an 
exclusive property of the Peierls transition. In fact, other 1-dim gapped spin models, as the 
XY chain with anisotropic interactions, present the same slowing down of the relaxation 
at special points (Ref. ^treats the Ising limit). A general connection between this critical 
behavior and the spectral properties is still lacking and currently under researchS. 

Our results are all exact, with closed analytical forms, from where the asymptotic regime 
for t — > cxD is obtained. This has been done using the fermion representation of the Jordan- 
Wigner approach for spin chains. Other procedures to solve the same systems are available, 
the most conspicuous being the Bethe Ansatz method, which applied to the XY model 
might hint at "hidden" conservation laws of the several versions of the model treated here. 
However, it is difficult to conceive a similar analytic calculation of the relaxation properties, 
using Bethe Ansatz techniques. 

We remark that integrable models displays other unusual dynamical phenomena, with 
anomalies in the transport propertiesil. For the Heisenberg model, the absence of spin 
diffusion has been probed by bosonization techniques^ and numerical calculation^!. Ref. 
m| presents exact numerical computations for the spin-1/2 XX Z chain at T = cxd, which 
include, as particular cases, the isotropic Heisenberg model and the XY regime, which are 
gapless, and the Heisenberg-Ising model, with a gap in the spectrum and long-range order 
for the ground state. Numerical data is inconclusive due to finite-size effects, but there are 
strong hints of a crossover from non-diffusive to diffusive behavior when we go from the 
gapless region to the gapped one. Concerning the low-temperature properties, a calculation 
by Sachdev and DamleE^ in the gapped region, shows that the diffusive behavior holds in 
the presence of long-range order. It is however puzzling to admit that low energy properties 
of the spectrum may affect the T = oo behavior. 

From a more fundamental point of view, we have studied the time evolution of a closed 
quantum many body system (at zero temperature), which is prepared in an arbitrary non- 
homogeneous initial state. The dynamics is solely given by the Schrodinger equation for the 
wave function or by the Liouville equation for the density operator (which includes the more 
general situation of a mixed ensemble for the initial state). Surprisingly, in the thermody- 
namic limit and for long times (A^ — > oo first than the limit t ^ oo), we get irreversibility 
in the form of a power law relaxation for the magnetization, in spite that the dynamics is 
unitary. No internal interactions are present, since all the Hamiltonians are reduced to the 
free-particle form. The irreversibility can be entirely ascribed to interference effects, that 
in our calculation are handled using the stationary phase method at asymptotic long times. 
Due to constructive interference, gapped models develop long-lived collective excitations 
with the texture of spin density waves (charge density waves for the fermion versions), that 
persist for longer times. We believe that those structures are universal features of gapped 
one-dimensional models. 

Finally, our simple quantum systems are good candidates for aging effects, a concept that 
has been coined to specifvnonequilibrium dynamics that depends on the initial conditions 
and relaxes very "slowly"Eil. But a complete characterization of quantum aging requires the 
calculation of the two-times correlation functionsS. 
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APPENDIX A: CALCULATION OF AVERAGES OF CANONICAL OPERATORS 



We illustrate our method with the calculation of the averages < ripf]^ >o, = ±) 

for models I and II. After taking the inverse of the Bogolyubov transformation 



< = Ukct + vtdt, fit = v^^ct - Ukdt , (Al) 



one obtains: 



< VpVg >o = UpUg < c'^c^ >o +v^v^ < d'^d'^ >o + 
+ UpV^ < c^di^ >o +v^Ug < d'^c'^ >o . 

Then using definition ( p!4D and relations 

where a" = + zaSj are the spin ladder operators and Lj = n/=i(25'J) is the string 
operator of the Jordan- Wigner transformation ([T^), it is possible to reduce the problem to 



calculation of the following spin averages: 

ia^„_i >o . (A2) 

All the averages above are taken with the initial density matrix p(0). If this is assumed to 
have the form (|llD, then most averages that involve products of different components of spin 
operators vanish. The only nonzero averages are those with products of 7— components of 
spin operators belonging to the same or neighboring lattice sites. Due to this property, all 
averages in ( |A2| ) vanish except for j = n,{n + 1). The results, for all the possible initial 
conditions p(0) = F{S^), are given below: 

a) case 7 = x 

< a2jL2jL2nO'2n >0 = < Cf2j-l-^2j-l-^^2n-l'32n-l >0= 2^j,n^fj,-a , 
< a2jL2jL2n-lO'2n-l >0= P^rtj < S2jS2j^i >0 ^ tt^nj+l < S2jS2j^i >0 ; 

b) case j = y: 



1 

— ( 
2 

< a2jL2jL2n~lO'2n-l >0 = —OiSnJ < 'S'ljS'lj-i >0 +Mnj+1 < 'S'fjS'lj+i >C 

c) case ■y = z: 

< a2jL2jL2nCt2n >0= ^j,nSfi,-a (^2~^^'^ ' 

< a2j_iL2j-lL2n-ia'2n-l >0= ^j,n5fj.~a (^2 ~'~ ^ '^2j-l ; 

< a2jL2jL2n~l<l2n-l >0= . 



< a2jL2jL2nO'2n >0 — < 0,2 j -1-^2 j~lL2n~lO'2n-l >0— -^^j,n^fj.-a i 
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Next, one calculates all the averages contained in the expressions for < Sqit) >o. We 
only give the results for 7 = x, ^ (the remaining case 7 = ?/ is identical with 7 = x): 

A) Case 7 = x : 

< VpVq >o= -{upv* + UqVp) < e'^iQ-) >o -[upV*e''' + UqVpC-'P] < >o , 

< PpPg >o= {upVg + Ugv;) < £^(Q-) >o +[upVge-'P + Ugv;e''^] < e^(Q+) >o, 

< ?7+/5g >o= {upUg - VpVg) < s'' [Q^ ) >Q +[upU gC'^^ " V pV gC" 'p] < s'^iQ^) >o , 

< PpVq >o= {upUg - v;v*g) < £^(Q~) >o +[upUge-'P - v;v*e''^] < £^(Q+) >o , 

where 

< e%Q^) >o= (< e^Q>o±< s^q^^ >o)/4M , 

< el >o= Ef=i < S^Sf^^ >o exp{tkj) ; 



B) Case 7 = 2;: 

< VpVq >o= lSp,g + UpUg < S'^iQ^) >o +VpV* < S^{Q-) >o , 

< p+pg >o= ^5p^g + v;vg < S'{Q+) >o +UpUg < S^iQ') >o 

< ri+Pg >o= UpVg < S'{Q+) >o -Mgt;p < S'{Q-) >o , 

< i3+r]g >o= Mpt;* < >o -UpVg < S'{Q-) >o , 

where 

< S^iQ^) >o= (< S^Q>o±< S^Q+. >o)/2M . 



Using expressions ( psf ) and (|24D , together with (|A3|) and (|A4|) , and taking the continuum 
limit when iV — *^ 00, we get the following expression for < S'q >4 : 

< S'q >t= ^ f {F^{k,Q)cos[tQt{k,Q)]+F-{k,Q)cos[m;{k,Q)]}dk, (A5) 

where ^ = 1, 2 , for the two models treated here. Quantities in (|A5|) are 

F^{k, Q) = ± 4:Ug{ul- I Vp \^)[t{e-'P/^ + e"^/^) < >o Im Vg 

for 7 = X , and 

F±(fc,Q) = A^{k,Q) < S^Q >o +5^(fc,Q) < 5^+^ >o , 

Q) = 1 ^ [1 - 2mJmJ - 2 I P -4MpMg Re {v;Vg)] , 
B^{k,Q) = ±AiUpUg Im {v*Vg) 

for 7 = 2r . 

Substituting the specific values (|17|) and (0) for {up, Vp) into formula ( [A5| ) one obtains 
the corresponding expressions for models I and II given in Section II. In the same way, 
analogous results may be obtained for model III. 
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APPENDIX B: AMPLITUDES AND PHASES 



We give in this Appendix, formulae for the amphtudes and phases of the oscillating 
components in the asymptotic development (|37D, to lowest order and for a generic Q-point : 



with 



and 



SiiQ,t) = a^,o{Q) [-j cos[tni{Q) + ai,o], (i = l,...,4), 



aifiiQ) = aloiQ) < S'q >o +Ko(Q) < Sq+^ >o, 
1 



«3,o(<3) 



"1,0 



1 



/Stt 
1 



57r 



(l + 5f(l-cosQ)/i?^+^, 



/2tt 
1 



= {l-6')smQ/Rl^^, g^±Qei 

TT 

(1 + cosQ), a^,o(Q) = 
(1 + cosQ), alJQ) = 



sinQ, Q^Qcs] 



sinQ, Qci <\ Q \< Q. 



c2 



j -n/A, 




1 ^/4, 





^2,0 







1 ^/4, 





(Bl) 



(B2) 



(B3) 



"3,0 



(7r/4)sign(Q) , a4fl = n/A 
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FIGURE CAPTIONS 



Fig. 1 Frequencies of tlie asymptotic components of the magnetization in (|37D , as 
functions of the wave number Q, for several values of the parameter 6. The different branches 
Qi, for i = 1,2, 3, 4, are indicated by the numbers. We also display the critical points iQci 
and ±Qc2 by solid squares. Note that the fl^ branch is tangent to the Qi and fl2 branches 
at and respectively, but it is only defined between critical points, Qd < |Q| < Qc2- 
The complete ^4 curve is shown as a guide for the eye. 

Fig. 2 Inverse of the relaxation times (tuiTj)"^ of the power-law decay of (|37D, as 
functions of the wave number Q for the same values of the parameter 6 displayed in the 
previous figure. The numbers label the different branches, in close correspondence with Fig. 
1. Solid squares are used to display (ci;ir4)~^, that vanishes at the critical points. The arrows 
in the upper figure indicate values out of the scale. 
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